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1. Introduction 

This paper is devoted to study some statistical aspects of the stabilization prob- 
lem for two and three-dimensional Oseen equations with help of feedback control 
defined on a part of the boundary that restricts a domain where the equations are 
determined. Since for stabilization problem the case of unstable equations is in- 
teresting, we assume that Oseen equations possess a solution that is exponentially 
growing as time t oo, i.e. the solution is unstable. New approach to the problem 
of stabilization by feedback control was proposed by one of the authors of this paper 
in [S|-mi- Namely, construction of stabilization from a part of boundary for para- 
bolic equation, 2D Oseen equations as well as for 2D and 3D Navicr-Stokes system 
was created in |H], 0, This construction reduces solution of a stabihzation 

problem to solving a mixed boundary value problem defined on an extended domain 
with zero Dirichlet boundary conditions and with initial condition belonging to a 
stable invariant manifold defined in an neighborhood of steady-state solution, near 
which we stabilize our system. In the case of (linear) Oseen equations, steady state 
solution equals to zero and invariant manifold is replaced on subpace invariant 
with respect to resolving semigroup and such that the solution going out initial 
condition wq e tends to zero as time i — > oo with the rate e""^* or faster. 
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Evidently, aforementioned mixed boundary value problem is not stable because 
if Wo ^ Xa then solution w{t, •) outgoing wq goes away from (and goes away 
from zero) even if wq is arbitrarily close to . That is why straightforward appli- 
cation of proposed construction of stabilization to numerical simulation may not be 
successful, because unpredictable fluctuations inevitably arise in real calculations. 
This situation was analysed in "TU", T3", and with help of the concept of real 
process introduced there. The method of damping of unpredictable fluctuations 
by some feedback mechanism was worked out in these papers and an estimate of 
stabilization for real process was obtained. 

This estimate is informative only when a norm of stabilized real process is not 
too small. But when this norm has the same order as norms of unpredictable fluctu- 
ations, aforementioned estimate became uninformative. Actually in this situation 
behavior of stabilized real process became chaotic. The goal of this paper is just 
to investigate the behavior of stabilized real process in small neighborhood of zero. 
More precisely, we solve the problem of retention of stabilized flow near unstable 
state of equilibrium. To do this we impose additional assumption on unpredictable 
fluctuations. We suppose that they are independent identically distributed (i.i.d.) 
random variables with probability distribution supported in a small neighborhood 
of origin for phase space. Then the real process is described by a random dynamical 
system and forms a Markov chain. 

Our aim is to prove that this Markov chain is ergodic, i.e. it possesses unique 
stationary measure invariant with respect to the corresponding Markov semigroup. 
This gives us the possibility to calculate by well known formulas the statistical 
characteristics of stabilized real process and to make clear its behavior using these 
formulas. 

During the last few years uniqueness of invariant measures for 2D Navier-Stokes 
equations have been proved by F.Flandoli, B.Maslovski [J], S.Kuksin, A.Shirikyan 
llll, ED, 123, W.E, J.Mattingly, Ya.Sinai jSj, S.Kuksin [TH|, Duan and Goldys |SI, 
and other authors. 

To prove ergodicity of indicated Markov chain we use recent results of S.Kuksin, 
A.Shirikyan jU] and S.Kuksin 18 on uniqueness of invariant measures for 2D 
Navier-Stokes equations with random kick-forces where coupling approach was ap- 
plied. 

Actually, it is enough for us to verify that random dynamical system arising in 
stabilization construction indicated above satisfies all conditions imposed in |21j . 
[18| on random dynamical systems. 

In section 2 we recall necessary information on stabilization method. In section 
3 we formulate the main results and present conditions imposed in j21| . jl8| on 
random dynamical systems in a form convenient for our situation. In sections 
4-6 we verify that these conditions fulfil for RDS corresponding to stabilization 
proceedure. 

We thank S.B.Kuksin and A.Shirikyan for useful discussion on ergodicity of 
Navier-Stokes equations. We thank also M.S.Agranovich for allowing us to read 
the proof of his result announced in ^ and formulated in Lemma l4.2l below before 
its publication. 

2. Preliminaries to the stabilization theory 

In this section we recall briefly results of |H1-C3 used below. 
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2.1. Formulation of the problem. Let ncR'^,d = 2,3, dn e C°° , Q = R+xU. 
We consider the Oseen equations: 

dtv{t,x) - Av + {a{x),\')v + {v,\')a + \'p{t,x) ^0, divv{t,x) = (2.1) 

with initial condition 

vit,x)\t=o=voix). (2.2) 

Here (t, x) = {t,xi, . . . , Xd) G Q, v{t, x) = (-ui, . . . , Vd) is a velocity of fluid flow, 
p(i, x) is pressure, a(x) = {ai{x), . . . , ad(a;)) is a given solenoidal vector field. 

We suppose that 9il = F U Fq, F 7^ where r,ro are open sets (in topology 
of dVl). Here, as usual, the over line means the closure of a set. We define S = 
M+ X F, Eo = X Fo, and set: 

t-lso-O, (2.3) 

where u is a control, supported on S. 

Let cr > be given. The problem of stabilization with rate a of a solution to 
problem (|2.1|I - H2.3|I is to construct a control u defined on S such that the solution 
v{t,x) of boundary value problem l|2.1(l - (|2.3(l satisfies: 

where c > depends on vq , a and Fq . 

2.2. The main idea of the stabilization method. Let w C M'' be a bounded 
domain such that D uj = ^, ^l^\uJ ^T. We set 

G = Int{nUuj) (2.5) 

(the notation Int A means, as always, the interior of the set A). 

We suppose that dG e C" where a > 2 is fixed and in all points except F\F = 9F 
it possesses the C°° smoothness. For the construction of uj and detailed description 
of dG in a neighborhood of 9F, please see (2], |13j . 

We extend problem H2.1|I " H2.2|) from il to G. Let us assume that 

a{x) eV'^iG)n {H^{G)y , a(a;) is real valued, (2.6) 

where, as usually, H'^{G), fc G N, is the Sobolev space of scalar functions, defined 
and square integrable on G together with all its derivatives up to order k and 
(i/''(G))'' is the analogous space of vector fields. Besides, i7g(G) = {f{x) G H^{G) : 
f{x)\x(^dG = 0} and 

V'^in) = {v{x) = {vu...,Vd)e (H^in))^ : div«-0} 

The extension of (|2.1|) - (|2.2|) from il to G can be written as follows: 

dtw{t,x) - Aw + {a{x),V)w + {w,V)a + Vp{t,x) = 0, divu;(i, a;) = 0, (2.7) 



w{t,x)\t=Q ^ wo{x), w\s = 0, (2.8) 

where S = M+ x dG. Note that, actually, wo from (|2.8(l will be some special 
extension of vo{x) from l|2.2|l such that the solution w{t,x) of problem H2.7|l . (|2.8|l 
satisfies the inequality 

\\w{t,-)\\vo^G)^ce-"'\\wo\\v«{G) fori^O (2.9) 
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For vector fields defined on G we denote by 751 tlie operator of restriction on f2 
and by 7r we denote tlie operator of restriction on F: 

7a : V\G) — > V\^), 7r : V\G) V^'-^'^iV), fc (2.10) 

Evidently, these operators are well-defined and bounded (see |23p. 

Definition 2.1. A control u{t,x) in (|2.1I) - (|2.3|) is called feedback ^ if 

v{t,-)=-inw{t,-), u{t,-)=-frw{t,-) (2.11) 

where {v{t, ■),u{t, •)) is the solution of stabilization problem (|2.1|) - (|2.3I) and w{t, •) 
is the solution of boundary value problem (|2.7() - (|2.8|l . 

Evidently, if the solution w of (|2.7() - (|2.8|l satisfies (|2.9|l . the pair {v, u) defined in 
(I2.11|l satisfies (|2.4|) . Since {v,u) satisfies (I2.1|l - (|2.3|l as well, it forms a solution of 
the initial stabilization problem (|2.1|I - H2.4|I . 

2.3. Description of "correct" initial conditions. First of all we describe the 
set of initial conditions {wq} such that solutions w{t,x) of (|2.7(l - (|2.8|l satisfy H2.9|l . 
Let G be domain 12.5|l and 

Vo'{G)^{vix)eV"iG):iv,,.)\an^O}, Vo^G) = V^G) n (H^oiG))' (2.12) 
where I'ix) is the vector-field of outer unit normals to dG. Evidently, 

MvO{G) = lkllv-«(G) — lkll(L2(G))<*; \Mv„HG} — ll'^^l! (L^CG))''^ ' 

Denote by 

n:{L,{G)r~^Vo'{G) (2.13) 
the operator of orthogonal projection. We consider the Oseen steady state operator 

Av = -ttAv + TT[{a{x),W)v + {v,V)a] -.V^iG) — >V^{G) (2.14) 
and its adjoint operator A*. These operators are closed and have the domain 
V{A) = V'^{G) n {H^{G))'^. Emphasize that V{A) consists of vector fields equal 
to zero on dG. The spectrums 'S{A),Y,{A*) of operators A and A* are discrete 
subsets of a complex plane C which belong to a sector symmetric with respect to 
M and containing In other words, A is a sectorial operator. So spectrums 

I](A), I](A*) contain only eigenvalues of A, A*, respectively. In virtue of (|2.t)|l they 
are symmetric with respect to M, and moreover S(A) = I](A*). 

We rewritten the boundary value problem H2.7() - (|2.8(l for Oseen equations in the 
following form 

^^^!^+Aw{t,-)=0, w\t=o = wo. (2.15) 

where A is the operator (|2.14ll . Then for each wq € Vq{G) the solution wit,-) 
of H2.15|l is defined by w{t, •) = e^^*wo where e""^* is the resolving semigroup of 
problem H2.15|l . 
Let (7 > satisfy: 

S(A) n {A G C : ReA = cr} = (2.16) 
The case when there are certain points of S(yl) which are in the left of the line 
{ReA = a} will be interesting for us. 

Denote by X+(A) the subspace of Vq{G) generated by all eigenfunctions and 
associated functions of operator A corresponding to all eigenvalues of A placed in the 

^It will be clear later why defined control really possesses feedback property 
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set {A e C : ReA < a}. By X^{A*) we denote analogous subspace corresponding to 
adjoint operator A*. We denote the orthogonal complement to X^{A*) in Vq{G) 
by X,{A) = X,: 

X„^VS{G)eX+{A*) (2.17) 

One can show that subspaces X'^{A), X^, are invariant with respect to the action 
of semigroup e"^*, and X„ + X+{A) = V^iG). 

Theorem 2.1. Suppose that A is operator (|2.14() and it > satisfies (|2.16() . Then 
for each wq € Xa- the inequality holds. Besides, the solution of problem \2.1f^) 

with such initial conditions are defined by the formula 



w{t, •) = e-^*wo = {2TTi)-' / {A - XI)-'e-^'wodX. (2.18) 

Here j is a contour belonging to p{A) := C \ S(^) such that argA = ±0 for 
A £ 7, |A| ^ N for certain 6 € (0,7r/2) and for sufficiently large N. Moreover, 
7 encloses from the left the part of the spectrum T,(A) placed right of the line 
{ReA = a}. The complementary part of the spectrum S(A) is placed left of the 
contour 7. 

Such contour 7 exists, of course. 
Proof: See 0, 

2.4. Theorem on extension. To complete the construction of stabilization for 
Oseen equations (12.11) . (|2.2(l we have to construct the operator E extending initial 
condition vq from H2.2(l from il in G such that Evq := wq G Xa-. This ^0(2;) we 
take as initial value in (|2.8|l . We consider here direct analog of construction from 

m 

Introduce the space 
V°{n, To) = Mx) e V°in) -. u ■ iy\r„ = 0, 3v e V^iG) : 70^} (2.19) 
supplied with the norm: 

ll"l|yo(n,ro) = inf \\w\\vO(g) 

w£V°(G):-inw=u " 

where 70 is the restriction operator defined in l|2.in(l and Vq{G) is defined in H2.12|l . 
We consider also the following closed subspace of Vq{G): 

V„°{G,n) = {weV„^{G): w\n=^Q}. 

For u £ V°{^, Tq) denote by u the element of the quotient space Vq{G) /Vq{G, 
such that for each v € u, jqv — u. Now we define the extension operator 

L: V^{n,Ta)-^V^{G) by Lit e u and Lu_LyO(G)Fo"(G, 17) (2.20) 

(i.e. Lu is orthogonal to Vq{G,^) in the space V§{G)). Evidently, ||i|| = 1. 

It is known (see jJOI) dl) that in the space X'^{A*) one can choose a 
basis {di(x), . . . , dK(x)) such that restriction {di{x)\uj, . . . , dK{x)\^) on an arbitrary 
subdomain u <Z G forms a linear independent set of vector fields. We can define 
space H2.17|l by the following equivalent form: 

Xa^{v{x)eVS{G) : ( v{x)-d,{x)dx = Q, j^l,...,K}. (2.21) 
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Theorem 2.2. f\l'2\. ) There exists a linear bounded extension operator 

E:V°{n,To)^X^, (2.22) 

i.e. (Ev){x) = v{x) for x G J7. 

Proof. Let subset lui C G \ il be a domain with C°°- boundary duji such that 
lnt{duji n dG) 7^ 0. In this set we consider the Stokes problem: 

—Aw{x) + Vp(a;) = v{x), divw{x) — 0, x £ toi] w\q^^ = 

As is well known, for each v € y°(ijji) there exists a unique solution w g V^^(tJi) n 
F^(a;i) of this problem. The resolving operator to this problem we denote as follows: 
(— 7rA)j^^i> = w. Extension of {~tt/S.)^^v from loi in G by zero we also denote as 
(-7rA)->. Evidently, (-7rA)-i^; € V^{G). 

We look for the extension operator E in the form 

K 



Ev{x) = {Lv){x) 



(x), (2.23) 



where L is the operator H2.20|l . cj are constants which should be determined. Evi- 
dently, Ev{x) — v{x) if X e for any Cj. Besides, Ev £ Vq{G). To define constants 
Cj we note that by (|2.21|) Ev £ X„ if 



dkix) 



K 

J2c,{-nA)zld,{x) 



dx^- dk{x){Lv){x))dx (2.24) 



for fc = 1, . . . , i^. As in [T2], ^3] one can prove that this system of linear equations 
has a unique solution. □ 

Aforementioned results imply the following result on stabilization (see [5]-|12)'). 

Theorem 2.3. Let domains f2 and G satisfy 1^2. 5\) . Then for each initial value 
vq{x) £ {fl, Tq) and for each a > there exists a feedback control u defined on E 
such that the solution v{t,x) of H2.1|I - (|2.4II satisfies the inequality 

ll«(ir)l!y„o(o) < ce-"* as t ^ (x. (2.25) 

2.5. Real processes. Recall results from 53 justification of numerical 
simulation of stabilization construction described above. In virtue of definition 
stabilization problem |jTIl)-(E3I) is reduced to problem (1^ . and ifT^ . 

Its numerical simulation is what we have to justify. 

Let e~^* = S{t) be resolving operator of problem (|2.7|l . H2.8|l . Suppose that we 
calculate this problem in discrete time instants 

ti <t2 < ...tk < ... 

where tk — kr and r > is fixed. Denote S — S{t). Let w'' be the result of our 
calculations at time instant tk. Since numerical calculations can not be exact, we 
have 

w'' = Sw''-'^ + / (2.26) 

where ip'' is an error of calculation which is unknown for us before time ifc. The 
sequence {w''} defined by H2.26(l is called real process. We suppose that we can 
estimate the error of our calculations a priori: 

y\\vOiG)^e«l, k>Q. (2.27) 



where £ > is a known quantity. Note also that at time tk the vector w'' is 
completely known (completely observable) since it is result of our calculations. 
Formulae H2.2(j|l is supplemented with initial condition 

w°=wo, woeX^. (2.28) 

(We assume for simplicity that <fi'^ = 0.) 
Equations ifT^ . ifT^ imply that 

k-l 

w'' = S'^wo + J2 SW^^' (2.29) 
i=o 

In virtue of H2.29|l . (|2.28l) the estimate ^ ce^^'^'^ is not true. Indeed, although 

S'^Wq G and therefore S'^Wq satisfies estimate of such kind, the fluctuations ip'^ 
possess nonzero components belonging to X'^{A). Hence, \\S^ •^'^'^WyOf^Q^ grows 
exponentially as fc ^ cx) that these terms destroy all stabilization construction 
described in previous subsection. 

To save this stabilization construction we use feedback mechanism (see |13 | .|14 | '): 
in the time moment tk when from H2.26I) is calculated we act on it by a spe- 
cial projection operator 11 : Vq{G) X^ that damps undesirable properties of 
fluctuations tp^ . This operator has to satisfy the following properties: 

i) Hip ^ip,yipe x„, a) \jif e vl\G) {Tiif){x) = ip{x) Vx e n. 

The operator satisfying these properties can be defined by the formula analogous 
to ^TT^ : 

y K -| 

Ilip{x) ^ ^{x) + Vc,(-7rA)jid, (x), (2.30) 



where constants Cj are defined by (|2.24|) with Lv — Lp. 

Applying to both parts of H2.26|l . H2.28|l operator 11. Taking into account that 
X„ is invariant with respect to action of operator S and using notation = , 
we get the recurrent formulae for controlled real process: 

lyfe+i ^ Sw^ + n¥j'''+\ for fc ^ 0; w° = (2.31) 

In ^21 1 d] the following estimate for controlled real process has been proved: 

Theorem 2.4. Suppose that unpredictable fluctuations ip^ satisfy \2.2T\l , and oper- 
ator n of projection on X„ is defined by 1^2.30^ . Then controlled real process {w^} 
defined by \2. 81]) satisfies the estimate 

ll^'llyo(G) < (7oe~""nko||yo(G) + ||n||£7(l - e"-^)) (2.32) 

where constant 70 G (0, 1) depends only on a if t > To and tq is a fixed magnitude 
^, ||n|| is the norm of operator ]2.'J(J\) . and i > is the magnitude from inequality 

(ITTtIi . 

Note that outside a small neighborhood of the origin, estimate (|2.32|l is equiv- 
alent to estimate (I2.9|l . Since in contrast to definition (|2.18|l of solution w{t, •) to 
problem l|2.7|) . H2.8|l . in definition H2.31|l of controlled real process fluctuations ip'^ 
permanently arise, therefore estimates (|2.32|l . H2.9|l can not be equivalent in a small 
neighborhood of the origin. Investigation of behavior for w'^ in this neighborhood 
is the main goal of this paper. 



'This property of constant 70 is obtained in formula 14.131 below. 
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3. Statistical problem in stabilization theory 



3.1. Primary considerations. Wc continue our investigation of the real process 
defined in subsection l2.5l above . 

It is clear from the definition l|2.26|l . (|2.28() of real process that its behavior is 
determined by its values in the instants tk — kr when the unpredictable fluctuations 
(fi'' arise. In the following we restrict ourselves by considering the real process only 
in these points and thus we obtain a discrete real process. So we consider iterated 
sequence (|2.31ll with a fixed initial vector wq = wo{x) G Vq{G) and unpredictable 
fiuctuation (p''+^ e V^{G) satifying fT?7|l . 

It is reasonable to assume that ip'' is a random variable (for each k) defined on 
a probability space {fl,A,F), taking values in Vq{G). We also assume that the 
random sequence (p'' is independently identically distributed (i.i.d.). Thus (|2.31|) 
defines a random dynamical system (RDS). 

For each fc, ip'' has distribution /i, whose probability measure is defined on 
the Borel ct— algebra B{V^{G)) of the space Vq{G), and is supported in a neigh- 
borhood of the origin: 

supp fi(lBi = {v€ Vo°{G) : \\v\\yo^a) < e}. (3.1) 

Let us consider the random dynamical system l|2.31(l from the point of view of 
stabilization. The solution of the stabilization problem can be derived in two stages: 
(a) To reach zero (or in general setting, to reach a steady state); (b) To keep the 
controlled solution w'' near zero. 

The solution of the stage (a) was explained above in Section 2. In particular, 
estimate (|2.32|l of stabilized real process in G was obtained. Here and hereafter, we 
denote — \\w''\\v^{G)^ unless otherwise noted. The stage (a) is the part of the 

process when w'^ does not reach the ball Brg with tq = |jn||e/(l — ce~°"^) where 
c e (0, 1) is the constant from H2.32|l . 

Note that if w'' ^ B,,„ , then using H4.7|l (see below) for estimating S we get 

WSw'' - w''\\ > \\w''\\ - WSw'^^W > (1 - 706--^^) II u;'^- II 

>(l-7oe— ) , ||n||g. (3.2) 

This estimate means that at stage (a), the distance between w'^ and Sw'^ is more 
than the norm of the fiuctuation Hip'', i.e., the deterministic component of random 
process w'' is prevailing. Therefore, in stage (a) the behavior of w''^^ is determined 
mainly by the term Sw''. 

The situation in stage (b) when w'^ e 5^ is different. Now the terms Sw'^ and 
Hep'' from the right hand side of H2.31|l have the equivalent order, and the motion 
of the realization w'' of the RDS H2.31|l becomes "irregular or chaotic" . 

The goal of this paper is to understand the behavior of the controlled stabilized 
sequence w'^ in the stage (b). We will do it with help of the modern theory of 
Markov chains or random dynamical systems. We begin with some preliminaries. 

3.2. Gauss measures. We recall some information about Gauss measures which 
will be used below. Let H he a Hilbert space with scalar product (•, •) and norm 
II • II, and let B{H) be the cr— algebra of Borel subsets of H. A measure G{du) 
defined on B{H) is called a Gauss measure if its Fourier transform G{v) is of the 
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form 



where a G H is the mean vector, and K : H ^ H is a hnear self-adjoint positive ^ 
trace class operator, called correlation operator of G: 

oo 

K* :^K>0, Ti&ce{K) = ^ < oo, (3.4) 

with {Xj} the set of eigenvalues of K. Differentiation of H3.3|l with respect to v 
implies that 



uG{du), {Kvi,V2) = {vi,a){v2,a) - J {vi,u){v2,u)G{du). 

Therefore a is the mathematical expectation of the Gauss measure G. 

In particular, when H ~ M™ and image ImK ~ M.™ , then for each T e S(R™), 

G(r) = / p{dy)dy, (3.5) 



where density p{y), y G M™, is defined by 

Piv) = (,,)J.detx - - (^-^^ 

Let A : iJ — > _ff be a continuous map in Hilbert space ff. As well- known ( 24 ), 
each measure /i induces a new measure A*fi{du) defined by 

A*fi{uj) = fi{A-^uj), Lu£B{H), (3.7) 

where A^^uj — {x £ H : Ax G lo}. This definition is equivalent to 

f{v)A*^l{dv) = J f{Au)^i{du) (3.8) 

for arbitrary / for which at least one side in equation H3.8|l is defined. Thus (|3.3|) 
and H3.8|) imply that if = G is the Gauss measure defined above and map A is 
linear, then A*G is a Gauss measure with mathematical expectation oi — Aa and 
correlation operator Ki — AKA*. 

3.3. Distribution of ip^ . We consider the right hand side of (|2.31|l where ip^ is 
an i.i.d. random sequence. We suppose that the distribution Vi^Lp^) of ip3^ has the 
form 

!?(/) = cxi{u)G{du) = v{du), (3.9) 
where c = (/g. G{du))~^ , and 

xe-(.)4;' ii^iif: (3.10) 

[0, ||u||>£ 

is the characteristic function of the ball defined in H3.1|l . G{du) is the Gauss 
measure with mathematical expectation a = and correlation operator K satisfying 
conditions H3.4|l 



■^Positivness condition can be weakened a little bit in the case of problem considered here: see 
below footnote in the subsection l5.2l 
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In virtue of definition H3.9|l . 

i^{uj) = cG{Bi Duj), Lje B{V^{G)) (3.11) 
and therefore z/(ijj) is a probability measure on B{Vq{G)) supported on the baU Bi. 

3.4. The main result. Since ip'' are i.i.d., Hip^ are i.i.d. as well, and D{Il(p'^) G 
Xa-- Thus RDS (|2.31|) defines a family of Markov chains in X^- with transition 
function 

Pik,wo,T)^¥{w''{wo)eT}, TeBiX,), (3.12) 

where w'' — 'w'^{wo) is defined by (|2.31|l and P is probability measure defined on 
(T-algebra A of subsets of probability space fl. Let VlX^r) be the space of Borel 
probability measures on B{Xc) and Cb{Xrj) be the space of continuous bounded 
functions on Xa-- Moreover, we denote by *P/c and the corresponding Markov 
semigroups acting in Cb{Xa) and ViXa), respectively: 

^kfiv)=^fiw\wo)) = J f{z)P{k,wo,dz), feCbiXa), 
^luiT)^ f P{k,wo,T)^l{dwo), ti^V{Xa), 

J H 

where E is for the mathematical expectation, and P(fc, wo,r) is defined in H3.12|l . 
A continuous function f{u) on X^ is called Lipschitz if 

\f{u + v)- f{u)\ ,. . ^ ^ 

sup — = lip/(u) < oo, M e Xa. 

■uex„,||D||<i 

We denote Lip/ = ||lip/(-)||c(x„)- 

A measure ji G V{Xa) is called a stationary measure for the RDS (|2.31|l if 
^P^,/i — fi,yk. The main theorem of this paper is as follows. 

Theorem 3.1. The random dynamical system \2.H1\) has a unique stationary mea- 
sure pi. Moreover, there exists a constant 7 G (0, 1) such that 

I j f{z)P(k, wo, dz) - J f{z)fi{dz)\ < C7^ A: = 1, 2, 3, • • • , (3.13) 

for every Lipschitz function f on X^ such that |l/|lc6(x„) ^ 1 '^i^d Lipf < 1. The 
constant c depends only on initial state ||it;o||- 

Note that the uniqueness of the stationary measure means that the RDS (|2.3HI 
is ergodic. The exponential convergence (|3.13|l means that RDS H2.31|l possesses 
the property of exponential mixing. 

Theorem l3.1l provides us the possibility of calculating easily the probability char- 
acteristics of Markov chain (|2.3HI . Indeed, in numerical simulation we actually 
obtain certain realizations w'' — w^{u}) of RDS H2.31|l . where a; G is a random 
sample. By the strong law of large numbers 

1 ^ f 
\\TCL ■j^'^'w^i^) ^ I wfiidw), asiV^oo, (3.14) 

where ix{dw) is an invariant measure of RDS (|2.3HI . In virtue of Theorem 3.1, 
li{dw) = jj,{dw). So while calculating w'^{uj), we can simultaneously obtain mathe- 
matical expectation of fi. Moreover, H3.13(l gives us the convergence rate in H3.14|l . 
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The topic connected with (|3.14|l will be studied in detail in some other place. 
Here we only note that strong law of large numbers was derived from ergodicity of 
Navier-Stokes equations in 19 in the case when random force is white noise. Of 
course in the case of kick forces this derivation can be made as well. 

3.5. Ergodic theorem. In order to prove Theorem 3.1, we use a result in |21[I18| 
on ergodicity of a Markov chain, proved by coupling techniques. Let us formulate 
this result in a form suitable for our problem. Consider a Markov chain (or a RDS) 
in a Hilbert space H 

u{k) ^ Tu{k ~ \) + rik, u{0) ^ uq, (3.15) 
where uq ^ H , T : H ^ H is & linear bounded operator such that 

llTull <7ohll, Vt^ei?, (3.16) 

for a constant 70 G (0, 1). 

Assume also that there exists an orthonormal basis {ej} in H (true for any 
separable Hilbert space H) and a sequence of subspaces Hi C H2 C ■ ■ ■ Hk C ... 
such that 

Hk = spanjei, . . . , Cr^} where — > 00 as k ^ co 

Denote by H^ the orthogonal complement to Hk in H : H ~ Hk ® H^, and 
designate by 

Qk ■ H — > Hk, Qk '■ H — > H^ 
orthogonal projectors. Suppose that 

\\QiTu\\<-fk\\u\\,yue H, and7fc ^ as A; ^ c» (3.17) 

At last assume that in l|TTK|l . rjk is an i.i.d. random sequence with distribution 
fj,(uj),uj € H, such that the projection Q^/i on Hk, has a continuous density R[x) 
with a compact support: 

Qlfi{dx) = R{x)dx 

where Hk 3 x ~ '^Xiet, dx = dxi ■ ■ -dxr^- Moreover, this density R{x) satisfies 
the condition: 



/ \R{x - vi) - R{x - ^2)1 dx c\\vi - V2\\Ht, 



(3.18) 



where the constant c > does not depend on vi,V2 G Hj^. 

Under aforementioned assumptions in | 21l I18| . the following theorem has been 
proved (see [TH|). 

Theorem 3.2. The RDS fff. _?5I) has a unique stationary measure fi. Moreover, 
there exists a constant 7 G (0, 1) such that 

I / fiz)Pik,u,dz)- f /(z)A(dz)| <C7^ k = l,2,--- , 



for every Lipschitz Junction f on H satisfying \f\ch(H) 5: 1 '"^'^ Lipf < 1. Here 
P[k,u,dz) is the transition function (see ^S.IS^) ) corresponding to RDS f3.15\) . 

4. Check of assumptions (|3.16I) . (|3.17() 

To prove Theorem 13.11 we have to check that RDS (|2.31|) satisfies conditions 
(|3.16|l . H3.17|l . and l|3.18|l of Theorem 13. 21 In this section we check the first and the 
second of them. 
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4.1. Subspaces of V^{G). First we introduce an orthogonal decomposition of 
Vq'(G) in order to define analogs of subspaces Hk from subsection 13.51 Recall 
that an important subspace of Vq{G) is (see definition in (|2.21|) ') 

X„ = {v eV^{G) : [ v{x)dj{x)dx = 0,j = !,■■■ ,m}, (4.1) 
Jg 

where {dj{x),j = !,■■■, m} is the basis of X^{A*), constructed in |^ from eigen- 
vectors and associated vectors of the operator A* , corresponding to all eigenvalues 
Xj satisfying ReXj > a. 

Let cr < (71 < • • • < (Tfc — > cxD as — > cx) be a sequence of numbers which satisfy, 
as a, the condition 

{XeC:ReX = at} D cr{A) =0 V fc. (4.2) 
Analogously to (|4.1(l . we can define the spaces 

Xa, - {v e Vo°(G) : / vix)djix)dx = 0,j = 1, - • • ,nfc}, (4.3) 

JG 

where basis {dj(x),j = 1, . . . ,71^} in X^^{A*) is constructed from eigenvectors and 
associated vectors of the operator A* , corresponding to all eigenvalues satisfying 
ReXj > ak- This basis is an extension of the basis in X^{A*) from (|4.1|) . and it is 
constructed by the same rules as the basis from (14. 1|) . 

Since ai > Uj > a for i > j we have that Ui > Uj > m and therefore X^. C 
Xa^ C Xa- Moreover, we introduce the following subspaces of Vq{G). Let X^ai 
be an orthogonal complement in X^ for the subspace X^-^, and Xg-j. 0-^+1 be an 
orthogonal complement in Xa-^ for the subspace ^crj,_,.i • In other words, Xa-a-i and 
Xa-kO-k+i are subspaces satisfying 

Xai ® Xaai = Xcr -'^fffc+i © ^CTkCrk + i = ^CTk- (4-4) 

We also define the subspace X^ C Vq{G), which is the orthogonal complement of 
X, in Vo^iG): 

X,®X^ ^V°iG). (4.5) 

Evidently, X^ — X^{A*). The subspace Xcr will play the role of space H in 
subsection 13. 51 Likewise, 

X,,^:=X„^,®---®X^^_^,, (4.6) 

will play the role of Hk and Xa-^ will play the role of . Recall that subspace X^, 
is invariant with respect to the operator S in RDS H2.31|l . That is why we put in 
(ITT5ll 

T = S\x,. 

Now we construct the basis {cj} from subsection 13.51 Let {ei, . . . , Cm} be orthog- 
onalisation of basis c?i, . . . , dm in X^{A*) = X^. Evidently, {ei, . . . , em} forms a 
orthonormal basis in X;^. Continuing orthogonalization process for 

{dm+l: • ■ ■ J dm },..., {c?nfc_i + l I • ■ ■ i '^Ufc }? ■ • ■ 

we get orthnormal basis {e,„+i, . . . , e„j} in X^ai and {e„fc_i+i, . . . , 6,^^} in Xa^:_-,ak 
for k — 2,3,.... We have to prove now that countable orthonormal system 
{ei, . . . ,ej, . . .} forms a basis in Vq(G). For this it is enough to establish that 
this system is dense in Vq{G). But in virtue of Keldysh Theorem (see J7|, |16|'). 
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the system {dj,j e N} constructed by eigenfunctions and associated functions of 
operator A* adjoint to operator (|2.14|) is dense in Vq{G). Hence, system {ej,j E N} 
obtained from {dj,j £ N} by orthogonahzation process is also dense. Therefore the 
system {ej,j = m + 1, . . . , m + fc, . . . } forms basis in the space H = 
In the next subsection we estabhsh inequahty H3.16|l . 



4.2. Certain properties of RDS ()2.31|) . We prove the following assertion. 

Lemma 4.1. Let S he the operator in \2.Sl]) and a > be given. Then for each 
7o G (0, 1), there exists r > such that for S = S{t) the following estimate holds 

\\Su\\<jo\\u\\,VueX„. (4.7) 

Proof. Recall that the basic space is Vq(G) and therefore we use the notation 
II ■ II = II ' llvo"- Inequality (|4.7|1 follows from the bound established in 

II / {A- Xiy^e-^-'dXW =\\ f {A + XI)-^e^^d\\\<ce-''^, (4.8) 

where A is the infinitesimal generator for S{t), and the contour —70- is defined as 
follows: 

where 

7^ = {A e C, ReX = -cr, ImX e [-{a + 0) tan(7r - -0), (cr + 0) tan(7r - ^p)]} 

7^ = {A € C, i?eA < -a, A = 7e±*'^ + 0, for 7 e [^^^,00)} (4.9) 

I cos tpl 

with 6 > and 7r/2 < -0 < tt fixed. We have to prove that c in (|4.8|l can be chosen 
independent of cr > 0. Since 70- belongs to the resolvent set of the operator A, we 
can get as in ,9| Lemma 4.7] that 

||(A/ + A)-i||<^^, AG -7., (4.10) 

with Ml > independent of A e — 7^- We see that 

iXI + A)-^e^^dX\\ s^h+h 

J --la 

where 

/i-/ ||(A/ + A)-i|||e^^rfA| h= \ ||(A/ + A)-i|||e^^rfA| (4.11) 
Jil hi 

Using H4.10|l we get 

<.((T+6l)tan(ir—i/') M.p-'^'^r^T i'(cr+0)ta.n{TT-ip) rl-r I n 



(cr+e) tan(;r-VO 1 + \/o^~fa? J - (ct+S) tan(7r-i/.) \/l + (x / a) 

.(l+e/fT)tan(7r-V>) j 

= Mie-"^ / ^ < Mice-"", 

where c does not depend on ct > 1 and r > 0. If ct G (0, 1), then 

(l+e)tan(^-V<) M^e-'^^dx 

J -{l+8)ta.n(7i-tp) ^ + Fl 
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with c not depending on cr e (0, 1) and r > 0. Moreover, by change of variables 

X = (7|cos-0| — 9), we get 

exp[— (7I cos^l — 9)T]dj 



h < 2Mi 



'■^Tj 1 + ^ (7I COS ipl — Oy + 7^ sin^ ip 



2Mi f°° exp[-a;T]da; ^ 2Mi f°° exp{~xT)dx 



IcosV'lJ,, 1 + ^/^^FTlxTl)HsIP^ \cosip\ 
2Mi f°° exp{-y)dy 2Mi exp(-crT) /"^ exp(-z)dz 



IcosV'lJtTT T + y Icost/'l Jo r(l + cr) + z' 

where in the final step, we used the change of variables y = z + ar. 
Note that 

exp(— z) dz 



T 



(l + cr) + Z 



1 if T $s 1 



For T S (0, 1) we have 



exp{—z)dz f°° exr>(—z)dz ^ dz , 

'-^ ' ^ ^ sC / h / e"^dz = l-lnr 



r(l + a) + z Jq t + z t + z 

So we have 

h < ^^e--^ci, (4.12) 
I cos V"! 

where ci does not depend on r > tq and a > 0. 

Thus we have proved that c in ()4.8|l does not depend on cr > and r > tq. So 

for u G Xcr, 

||5'(r)u|| < cAfie"'^^. (4.13) 

where c > does not depend on ct > and t > tq. Therefore for given 70 and 
cr > 0, we can take r > such that ce^°"^ = 7o- O 

Inequality H3.16|l evidently follows from H4.7() . 

4.3. Estimate S{t) on Xcr^ for large k. We check here bound (prT7|) for T = S'(t) 
and Q^H — Xa-^. Simultaneously we make more precise the choice of ak- Together 
with operator A defined in (|2.14|l we consider the Stokes operator 

Ao - -^A : Vo°{G) ^ Vo°{G) (4.14) 

where ir is projector (|2.13|) . Operator Aq is positive self-adjoint with domain 
T>{Ao) = ViA) = V'^{G) n Vq{G) and with discrete spectrum. Let {ej} be eigen- 
vectors of ^0 forming an orthonormal basis in Vq{G) and < /ii ^ /i2 ^ . . . 
be corresponding eigenvalues, taking into account of their multiplicities. Then for 
each g e M the power A^ can be defined by the formula A^v = ^ ^j)v°(G)^j- 
Eigenvalues nj possess the following asymptotic for large j: 

H^Pof' + Oif'^llnj) as j ^ 00 (4.15) 

where /3o > 0, d is dimension of G (i.d. d = 2 or 3), 0{j'^^'^/hij) is a function sat- 
isfying the estimate \0{p^'^/hij)\ ^ cj^/''/ In j as j ^ 00 with c > independent 
on j. Asymptotics (|4.15|l was obtained by K.I.Babenko 0]. 



^Actually, K.I.Babenko proved 14.151 in 4 in the case d=dimG=3 only. But his proof can be 
extended in the case d=2 as well. 
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We can write operator A from (|2.14fl as follows: 

A = Ao+Ai (4.16) 

where Aiv — Tt[{a{x),W)v + (w,V)a]. Using |15[ Ch.3, Lemma 4.5] one can easily 
get the following bound: 

\\AiA-'/^\\ =b<oo (4.17) 

Now we describe one result of M.S.Agranovich announced in Bound (6.61)]. 
Although it is obtained for general abstract operators we formulate it in the case of 
Oseen operator. Its proof will be published in 2 .This result consists of the choice 
of sequence ct/c — > cxd as fc — > cxd such that on segments 

Tfe = {A = (Tfe + ij; I7I ^ b'al^^}, b' > does not depend on fc, (4.18) 

resolvent {A — A/)^^ possesses some optimal bound. Numbers ak are found on 
segments ^ 

Afc = [e^'^/^ e^e^+i)/''], d = dimG 2 or 3. (4.19) 

Lemma 4.2. Suppose that an operator A has the form i4.1b\ l where Aq is self- 
adjoint positive operator with discrete spectrum and eigenvalues satisfying i4-.15{ l, 
and Ai satisfies j^.j7| ). Then there exists Uk G such that for A e Ffc the 
following estimate holds: 

d-l 

\\{A^\I)-^\\^c,c''/ al^l^ (4.20) 

where segments F^,, are defined in \4.1^ , \j^.iy'^ , respectively. 

Using Lemma r4.2l we can check (|3.17ll . Recall that S[t) = e^"** is the resolving 
semigroup of problem (|2.15|) where A is Oseen operator (|2.14ll . and S = S{t) where 
r is a fixed number chosen in Lemma [4 . 1 1 such that H4.7|l is true. Note that if we 
would increase r, H4.7|l is true as well. Recall that the space Xa-^. defined in ()4.3|l is 
invariant with respect to the operator S = S{t). 

Theorem 4.1. Let A be Oseen operator S — S{t) — e~^'^ , and sequence 

(Tfc 00 as k 00 be chosen in Lemma \4.^ Then there exists Tq > such that 
for each t > tq on the spaces the following estimate hold: 

\\S{t)u\\x,^ ^lk\\u\\x„^^, where 7fe as k ^ 00, (4-21) 

where 7^s do not depend on u ^ ^a^- 

Proof. It is clear that for u e X^r^ 

S{T)u^{2m)-^ I {A-\I)-^ue-^''d\ = -{2m)-^ I [A + \iy 
where —70-4, — 7^^, U 7^^, and similarly to 1)4.9(1 

7^^ = {A e C,i?eA = -ak,Im\ e [-(cr^ + 61) tan(7r - (cr^ + 61) tan(7r - -0)]} 

7^, - {A e C, i?eA < -a,, A = 76^^'^ + for 7 e [f^, 00)} 

I cos?/;! 



ue 



^^d\ 



^Actually, Agronovitch's result was proved in |5| under assumption that reminder term in 
I4.15i has the form 0{j^) with r < 2/d (that is stronger than 0(j^''''/ In j) from i4.15i 'l. and 
are looked for in segments [k'^P/'^, [k + l)2p/''] with some p > 4. But if we change these segments 
on segments 14.191 then a streightforward repeating of the corresponding proof from 1^ leads to 
the desired estimate. 
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Doing calculation as in proof of Lemma 14.11 we get the same formulas where cr is 
changed to Ufc only. Then estimation of the term I2 gives as in ()4.12f) : 

h < I lie-'^' ci, 4.22 

I cos V'l 

where ci = ci{tq) does not depend on r > tq. To estimate Ii we use (I4.20|) instead 
of (|4.1U|) . More precisely we use the estimate 

||(A + A/)-iKcie"^''"^V;'/', Xe-fk (4.23) 

where ffc = {A = cr/j + ij; \j\ < {ak + 9) tan(7r - ^)}. If A G -Ffe, estimate (I4.23f) 
directly follows from H4.20|l . For A e — {F/c \ F^.} situation is easier and one can 
prove estimate similar to (|4.10|l in right side of which |A| is changed on |A|^/^ that 
is stronger than H4.23|l . (This has been done in I2|). Applying H4.23|l to I2 defined 
in (|4.11|l (with a changed on ak) we get: 

/((Tfc+fl) tan (tt— 
-(o-fc+e)tan (tt-i/)) 

where c does not depend on fc. If we choose t > lnc2 then 14.22|l . (|4.24|) imply that 

h + h ^ ^ as k ^ CO. 
This proves (jOU- □ 
Inequality H3.17|l follows from (|4.21|) . 

5. Reduction to the finite dimensional case 

Now we need only to check the condition (|3.18|) . The rest part of the paper is 
devoted to prove that the distribution Hip'' of random forces in RDS (|2.31() satisfies 
this property. This will then complete this paper. First, we project Il(p'^ on finite- 
dimensional subspacc. 

5.1. Calculation of the projection for probability distribution Hip''. In this 
subsection we calculate probability distribution for random variable Sw'^ + Hip'^'^^ 
in RDS (|2.31|) under assumption that Sw'^ is a fixed vector. Since Sw'' € Xcr and 
n : Vg" — ^ Xa- is a projection on the probability distribution 'D{S'w^ + Hip^'^^) 
is supported on X^- It is enough to calculate Vlliip^^^) because ^{Sw^ + liip^^^) 
is simply the shift of Vijiip^^'^) along the vector Sw^ . In virtue of H3.9|l - H3.11|l . 
T>{tf*'^^) = where the measure v is defined by 

v{ijj) = cG{Bi n w), e B{VS), (5.1) 

where G is the Gauss measure with mathematical expectation a = and correlation 
operator K satisfying (|S31l. Clearly, V{TlLp^+^) = W V. Since 

H:V^ ^ X, (5.2) 

is a linear projection on Xg., we have by (|3.8|l for uj E B{Xrr): 

U*i^{uj) = c /" x^(nu)G'(du) = c /" G{du) 



c / G{du) = c / n*G(dw). (5.3) 
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We have already mentioned that II* G is the Gauss measure supported on Xa- with 
mathematical expectation a = and correlation operator Ki = HKH* . Note that 
in (|5.3|l . Xi^(s') is the characteristic function of the set u, and 

n-^uj = {xe : Hx e u}, fi-^uji = {x e n Bg -. nx e t^i}. (5.4) 

Hence to define completely the measure li* v{ijj) from 1)5.3(1 . we need to calculate 
liBg. Let us consider the decomposition 14.5(1 of Vq*. Since 11 : Vq* ^ C Vq* 
is a linear projection, for each y € Vq decomposed as y = y' + y" with y' E X^-, 
y" E X^, we obtain 

X ^Uy = y' + Ay" where A = Il\x±. (5.5) 

So 

A : X^ Xcr, A* : X^, — > X^, 

where A* is the operator adjoint to A. (We identify Hilbert spaces Xcr,X^ with 
their dual spaces.) Note that 

B. = {\\yr<i'} = {\\yT + \\y"r<e'}. 

Thus by CTl . 

TlBg = {xeX„: There exists y" E Bg n Xj-such that \\x - Ay"f + \\y"f < e^}.{5.6) 
To make this more precise, we consider the extreme problem 

f{y")^\\x-Ay"r + \\y'T^ini, y" E X^ . (5.7) 
The solution ij of this problem exists, is unique and satisfies 

(/'(y), /i) = 2{-(.T - Ay, Ah) + (y, h)} - 2{y + A* Ay ~ A*x, h) = 0, Vh E X^. 
Since the operator A* A is nonnegative, the equation 

A* Ay" + y" = A*x 

for each x E X^r has a unique solution y — y{x) E X^ and it is the solution to the 
extreme problem ((5.7(1 . It is clear that the map x ^ y{x) : X^ X^ is a bounded 
linear operator and y{x) — [A* A + E)~^A*x. Thus the definition of XlBg in ((5.6(1 
can be rewritten as follows 

nSg ^{xeX„:\\x~ A*y{xW + \\v{x)f < e'}- (5.8) 

The set ((5.8(1 is an ellipsoid. Thus 1(5.3(1 and 1(5.8(1 define the measure n*i^, and 
hence the RDS ((2.31(1 is completely defined as well. 

5.2. Finite-dimensional measure. In order to prove that the RDS 1(2.31(1 is er- 

godic, we study a map of the measure 11* i^. We introduce the operator of orthogonal 
projection Q connected with projection operator Qk from Subsection 13. 51 

Q:V°-^X^^~X^®X^,„ (5.9) 

where subspaces X;^, X^ak are defined in ((4. 5(1 . ((4. 6(1 . We have the following prop- 
erty for the operator Q. 

Lemma 5.1. Let 11,(5 projection operators defined in 1^2. 3 U\) and \5.y\) . respec- 
tively. Then 

gn = gng. (5.10) 
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Proof. Each x G Vq can be decomposed as follows 

a; = yi + 2/2 + 2/3 

with J/1 e Xj-, 2/2 e -'i^CTCTfc and j/3 e X^-^ . Then 

Qnx = Q{nyi + 2/2 + ys) = Q^yi + y2, QUQx = Q{nyi + ys) = Qn^i + ys- 

□ 

The random dynamical system (|2.31() generates naturally the measure n*i/. We 
now study the measure Q*Il*i/. By H5.10|l we show that this study can be reduced 
to the study of a measure defined on a finite dimensional space. 

Theorem 5.1. Let Q be the orthogonal projection in 11 be the projection in 

\2.S(]\) . and v{ijj) be a probability measure on B{Vq). Then 

Q*U*iy ^Q*U*Q*i^. (5.11) 
Proof. By (ESJ and we get 



f{u){Q*n*iy){du) ^ J f{Qnv)i^{dv) 

fiQUQw^dw) = J f{u){Q*U*Q*iy){du) (5.12) 
This proves the theorem. □ 

The relation H5.11|l allows us to reduce our investigation to the case of measures 
defined on finite dimensional space. Let us calculate the measure Q*!^. Taking into 
account the definitions (|3.9|) of ^{du) and (|5.9|l of Q, we get, analogous to H5.3|l . 
that for each T e BiX^J, 

Q*v{r) = c [ Q*G{du), (5.13) 

where Q*G is the Gauss measure supported on X^^, with mathematical expec- 
tation zero and correlation operator QKQ. We make the following identifica- 
tion taking into account identity in (|5.9|) : using in X^^ an orthonormal basis 
{cj , j — 1, . . . , n} (n = rik) introduced in Subsection 14.11 we see can write 

n ra n 

and take the following identifications: 

X^^ - M" = {y = (yi, . . . ,y„)}, - R" = {u}, X,,, - M"-" - {v} . (5.15) 

We restrict correlation operator QKQ on X^^. Then QKQ can be regarded as a 
n X n matrix. By (|3.4|l this matrix is non-degenerate because for each 7^ u G 
{u,QKQu) = [Qu.KQu) = {u,Ku) > and therefore kerQKQ = ^ We 
denote K = (QKQ)'^. By (E^J-lESJl, we conclude that 

Q*G{dy):^G{dy)^g{y)dy, (5.16) 



^Actually we can weaken the first condition in 13.41 assuming that K* = K ^ and ker QKQ - 
where Q is orthogonal projector on ^i^^, with big enough ak 
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where g(y) = ii£t4e-^(-^>'2'). Note that 

n n 

We hence obtain 

P{dy) := QV(dy) = cXe(y)G(dy), (5.18) 

where c = (/^ G{dx))~^ and Xe is the characteristic function of the ball B. 
We introduce the projection operator 

TT = QU : X^^ ^ X,,„ (5.19) 

In virtue of (|5.18|) . for each oj e B{X^J, 

g*n*QV(w) = 7r*i>(w). (5.20) 

6. Density P{x) of the measure 7r*j> 

6.1. Preliminaries. Recall that the space X^^ admits the orthogonal decomposi- 
tion 

X^^=X^®X,,,. (6.1) 

Below in order to emphasize belonging u E X^ , v € ^o-o-j. we write u(Bv instead of 
u + V. 

The projection operator tt defined in (|5.19|) can be represented as follows 

7r^ia,E), (6.2) 

where E is the identity operator in X^-a-^ and a : X^ ^a-ak ■ For x G we 
denote by tTx , the afhne plane in X^^ : 

TT^ = TT^^x = {y = u®veX^® X„„^ = Xj-^ ■.au + v = x}. (6.3) 

In particular, when a; = 0, 

TTo = TT^^O = {2/ = it®weXj^e X„„^ ■Ciu + v = 0}. (6.4) 

Since B is the support of the measure v{dy) = v{du,dv) defined in H5.18|l . the 
ellipsoid ttB is the support of the measure 7r*z). 
For each / G C{ttB) we have 

f{x)TT*i){dx)— / f{au + v)C'{du,dv)— / f{x)dx / r(u', a;)du', (6.5) 

B Jb JttB J-n^nB 

where the first equality is via the definition of the measure Tr*^' and the second 
equality follows from the change of variables 

w ^ u — au G ttq, X = au + v (z X^^. (6.6) 

The calculation of r(?«, x) will be done later. The formulae (|6.5|) gives the expression 
of the density P{x) for the measure T:*i>{dx): 

P{x)dx — TT*i'{dx), where P{x) = / r{w,x)dw. (6.7) 

JTT^nB 
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6.2. Change of variables {w,x) — > {u,v). In order to calculate the kernel func- 
tional T(w,x) in (|6.7|) . we need to consider the following change of variables, i.e., 
the inverse of (|6.t)|l 

u = u(w,x), v = v{w,x). (6.8) 

We introduce an orthonormal basis {bj, j = 1, • • • , n} in X^^ . Let 

{bj,j — 1, - ■ ■ , m} be the orthonormal basis of C X^^ , (6.9) 

composed of eigenvectors of the operator E + a*a : X^ X^. Suppose that 1 < 
/^i, ■ ■ ■ : 1 ^ A'm are eigenvalues corresponding to eigenvectors 6i, • • • , bm- Assume 
that 

Ml > I7 ■ • ■ , Ms > 1, Ms+i = • • • = Mm = 1- (6.10) 

Lemma 6.1. The following statements hold: 

(i) bj e ker a, j = s + 1, • ■ ■ , m; 

(ii) {abj, J = 1, • • • , s} form an orthogonal basis in Ima C X^^^. 

Proof, (i) For j s+1, • • • ,r7i, {E+a*a)bj ~ bj iSa*abj ~ 0. Since ker a* _L Ima, 
abj ^ implies a* abj 7^ 0. Hence bj G ker a. 

(ii) Note that for i,j = 1, • • • ,s, (a6i,a6j) = {a*abi,bj) = {fii — l){bi,bj) = 
and \ii ^ \. Hence abi _L abj. If 7 G Ima then for a certain 6 G ^i", we have 
7 = a6 = "X^jli "^j^j — Si=i ^i'^^i- Therefore {a5j, j = 1, • • • , s} form a basis in 
Ima. □ 

Since ker a* Ima — Xa-a-f. and dim Ima — s due to Lemma 16.11 we see that 
dim ker a* — n ~ m — s. Let ■ • ■ , b„} be an orthonormal basis for ker a*. 

Then by Lemma l6. II again, the vectors 

ab ' 

^ TTT^'-^ = 1' ■ ■ ■ bm+s+i,- ■ ■ ,bn (6.11) 

\\abj\\ 

form an orthonormal basis in Xcrak- In virtue of (|6.9|I - H6.11|I . 

Vectors 61, • • • ,6„ form an orthonormal basis of X^^. (6-12) 
On the plane ttq defined in H6.4|) . we consider the vectors 

{^ + \\abj\\^)^ 

Lemma 6.2. Vectors \6.13\) form an orthonormal basis of the plane ttq. 
Proof. By the definitions H6.2II . (|6.4|l and (|6.13|) . we see that 9j G ttq, as 



ne,^{l + \\ab,r)-n^b,^ab,]^0. 



By Lemma |6. II for i =/= j, 



(ft a\ - ,h®{-ab^ bj®{-abj) 

= (1 + + ||a6,|p)-^[(6„6,)^^^ + -a6,)x„.J 

= (6.14) 

and therefore the system (|6.13|l is orthonormal. As the rank of the matrix for (|6.2|) 
equals to n — m, 

diniTTo = dim ker TT = dimker(a, E) — m. 
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Hence the system H6.13|l forms an orthonormal basis of the plane ttq. □ 
Define 

9j — bj, j — m + I, ■ ■ ■ ,n. (6.15) 

Then the vectors 

0„ J = 1, (6.16) 

defined in (|6.13() and (|6.15() form a basis of X;^^ . 

Let R ~ (Rij) be the n x n matrix with components Rij defined as follows. 

Ru = (1 + \\ab,\f y^; = , " , , = , i = 1, • • • ,s; 

VI + \\ah\\^ 

Rii = 1, i = s + 1, • • • ,n; Rij = for other (6.17) 

It can be checked directly that the matrix R transforms the basis b — (fei, • • • , 6„) 
to the basis 6 — {di, ■ ■ ■ , On)'- 

9 = Rb. (6.18) 

Now we can calculate the change of variables H6.8fl . the inverse of (|6.6f) . Let 
y € X;^^ admits decompositions y = u + v = w + x with u G X^ , v G ^o-o-^ , w G ttq 
and a; € ^o-o-fc • Define 

n n 

i=i j=i 

m n 

j—l j—ra+l 
m n 
W — ''^^WjOj, X= Xj-mdj- (6.19) 

j—l j—7n+l 

We introduce notations 

y = {yi,--- ,yn),z = {zi,--- ,z„), 

u = {uj =yj,j = 1, - ■ ■ ,m),v = {vj = yj+,n,j = 1, • ' ' , - m), 

w ^ {wj = Zj,j = 1, • • • ,m),f = {xj = Zj+„iJ = 1, • • • ,n-m). (6.20) 
Then the change of variables H6.8|l is rewritten as 

u = u(w,x), v — v{w,x), (6.21) 

or in more compact form 

y = y(z). (6.22) 

Theorem 6.1. The transformation j6'. 21^) can be calculated as follows 

y^y{z)^R*z, (6.23) 

where R* is the conjugate matrix of the matrix R defined in Note 
that the matrix R transforms the basis {bj} to the basis {9j}. 

Proof In fact, the relation (f^:^ follows from □ 
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Using (|(j.23|) and (j6.17[l . we obtain the Jacobian of the transformation (|6.23|1 : 

s 

J = detDy/Dz = detR* = + Hab^f )-|. (6.24) 

i=l 

We see that the Jacobian J depends only on the operator a. Now we make more 
precise the expression for density P{x) in (|6.7|) . Note that in l().5|l we made just 
the change of variables l().22|l or the equivalent H().23|l . Taking into account of the 
definition y ~ {u,v) and z — {w,x), and using the facts (|5.16() . H5.18|l . (|6.23l) and 
(|6.24|l . we obtain that the integrand T{w,x) in the expression of density P{x) in 

s 1 

T{w,x) = YI{1 + \\abif)i detk{27r)-"/^ exp[-~-{kQR*z,QR*z)], (6.25) 
with z = {w,x). We suppress the arrow — > on top of w and x here. 

7. Smoothness of the density P{x) 

Formulas (|6.5I) and H6.7|l imply that the density P{x) is supported in the ellipsoid 
ttB C Xaa^ ■ So P{x) = for o; G X^a^, \ T^B. We now investigate the smoothness 
of P{x) for X £ d{TiB) and for x E Int(7ri?), respectively. 

7.1. Smoothness of the density P{x) on boundary d{TTB). It is clear that the 
set n TTj; is if X ^ ttB, it is a single point if a; G d(nB), and it is a ball in the 
n — ra dimensional plane tTj^ if a; G Int(7ri3). We first calculate the center and radius 
of this ball. 

Lemma 7.1. Let the hall B be defined in \5.n\j , x G Int^irB), and the plane tt^ be 
defined in <6'.5^) . Then the center of the ball B H tTx is 

w = u®v = [{E + a*a)~^a*x] ® [x - a{E + a*ay^a*x] (7.1) 

and the radius of the ball B H tTx is 

r = (e^ - \\{E + a*a)-'^a*xfx^ - \\x - a{E + a*a)-^a*xfx^^^)K (7.2) 

where, recall that, e is the radius of the ball B. 

Proof. Evidently, the center {u, v} is the solution of the extreme problem 

\Mx^ + IkllL., ^ iiif' ^ (7-3) 

By definition H6.3|) . {u,v} G tTx iff au + v = x, i.e., v = x — an. Substituting 
this into (|7.3|l and solving the extreme problem, we obtain the solution (|7.1|l . The 
radius r follows from the Pythagoras theorem. □ 

Let us consider the following extreme problem: Given x G ttB, find h G X„„^ 
such that 

F{h) = \\{E + a*a)-^a*{x + K)f + \\{E ^ a{E + a*a)-^a*){x + h)f inf, 

(7.4) 

\\h\? = ll (7.5) 
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with 7o > a given sufficiently small parameter. Recall that each x G ^o-o-fc admits 
the decomposition 

s 

X = xq + xjabj, xq G kera*, xj e M, (7.6) 

where {bj} is the basis in Ht).9|l . Note also that, by Lemma Iti.lf ii). {abj,j = 
1, • • • , s} is a basis of Ima. 

Lemma 7.2. Suppose that x G ttB and it has decomposition ^7.6]) . // 70 > is 

small enough, then there exists a unique solution h of the extreme problem i7.4{ l- 
i7.S}) . The solution h is determined by 

s 

h = ho + hjabj, ho G kera*, h^ G M, (7.7) 

3 = 1 

where 

l + /(7o) 1 + /(7o)Mj 

and xo,Xj,j — 1, . . . , s are defined in i7.ti{ l, fij > 1, j — I, ■ ■ ■ , s are eigenvalues 
f6.10\} of the operator E + a*a, and ^(70) is the unique solution of the equation 

{i + iY^j-^{i + it.,Y-^'- ^'-'^ 

Proof. The existence of a sohition of the finite-dimensional problem H7.4(l - ()7.5|l is 
evident. We now prove the uniqueness of this solution h. 

Let £{h,l) = F(h) + /(||/i|P — 7o) be the Lagrange function for the extreme 
problem (|7.4|) - (|7.5|) . By the Lagrange principle, if /i is a solution of this problem, 
then there exists I G R such that 

{C'h{h,l),5) = {F'{h),5) + 2l{h,5) = 0, V5 G (7.10) 

Substitution of expression F{h) in (|7.4I) into (|7.1UI) yields 

{{E + a*a)~^a*{x + h), {E + a*a)-'^a*5) 

+ {{E - a{E + a*a)-^a*){x + h), {E ~ a{E + a*a)-^a*)5) + l{h, 5) = 0. 

We transform operators from right multiplies in scalar products to the left multi- 
pliers. Noting that 

a{E -I- a*a)-'^a* + {E - a{E + a*a)-^a*f ^ E - a{E + a*ay'^a*, (7.11) 

we obtain equations for h and I: 

{E - a{E + a*ay^a*){x + h) + lh = 0. (7.12) 

Substituting the decompositions (|7.6|) - (|7.7|l for x and h into (|7.12l) . we arrive at the 
following system of equations 

xo + ho{l + I) = 0, (7.13) 



[1-^ -]{xj+h,) + lh,=0, 
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(7.14) 



and therefore 

By (I731),l(77|) and lf7?T5|l . we finally get 

- ii/^ip = whor + E /^iii^fe.ip = ^ + E 1^ = ii (7.16) 

Since r(0) = ||a;||^, r{l) ^ as / ^ oo and r'(Z) < 0, there exists a unique solution 
/-;(7o)if7o'<ll^ll'- □ 

Theorem 7.1. Let x e ^ttB, x + he Int{nB) and \\h\\ < \\x\\. Then P{x) = and 
for some positive constant c, 

\P{x + h)\<c\\h\\'^ as\\h\\->Q. (7.f7) 
Proof. If X e di-KB), then 

where the point u(Bv is defined in (|7.I|) . Hence by (|6.7|l . P(a;) = 0. Denote = 70 
and take the solution h of the problem im|) - (|731) . By (fT^ . lfrU) and lf73|l . the 

radius r{x + /i) of the ball B fl tt^^/j is maximal in the sets of radius r{x + h) of 
-B n 7ra;+;i corresponding to vectors h such that \\h\\ = 70. We calculate r{x + /i). 

Since x £ d{TTB), we have ||u(a;)|p + ||'0(x)|p — with {u,v} defined in H7.I|I . 
In H7.2|l . taking x ^ x + h and substituting = + ||-C(a;)|l^, and taking 

into account of H7.6() - (|7.8() . we get 

r^{x + h) = -{{E + a*ay^a*h, (E + a*ay^a*{2x + h)) 

- {{E - a{E + a*a)-^a*)h, (E - a{E + a*a)-^a*){2x + h)) 



{{E + a*a)-'a* E (E + a*ar'a* E ^i^H^ ^ — ^)) 

j = l r J 

(2 - ^) + E(^^ - + a*a)-^a*)a6zxi(2 - j^)) 



- 1) xi{fj.i - 1) ^^ _ _J_^(h h \ ^ ll^"lr(^ + ^0 
^u,(i + 1 + ;^/^'^'''^+ (i + o^ 



+ E xMab,,ak){l - ^^^)t-^(1 - )(2 



_^ ^j(Mj"l)'(l + 2;/i,-) ||a:o|p(l + 2/) ^ 2 , ||2 1 + 2?M, 

(7.18) 

We rewrite this as follows: 

, . ted^ , x: ,f ' " (1 + (7.19) 
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There exist constants < ci < c^/^ such that for each j = 1, • • • , s and for every 
Z > 0, we have 



ci(l + < + ^^Vrf) < ^(1 + 0- (7.20) 



Comparing H7.9|l and H7.19|) . we thus get 

ci7o(l + < r''{x + h)< C27o(l + 0- (7-21) 
Let Ai{x) = (ll^olP + ELi Then using ^ we get 

and therefore 

l + /(7o)>^. (7.22) 

70 

Let now ^2(0;) = (Ijxof + Ej=i ^U^^J H^)^- ^^^^ obtain 



7o 



2 < ^2(2;)2 



Hence, 

l + ^(7o)<^. (7.23) 
70 

Substituting |(021l, into ((ZiHl), we get 

cMxho <r^x + h)< C2A2{x)jo- (7.24) 
By (|6.25|l . there exists constants < ci < £2 such that 

ci <T{w,x) <C2, (7.25) 
for each {w,x) — z such that ||Qi?*z|p = ||t/|p < e^, i.e., on the ball B. Hence, by 



ciV{Tr^+hnB) < P{x + h) ^ T{w,x)dw <C2ViTT^+hnB), (7.26) 

where ViiTx+h H B) is the volume of the ball tTx+ii H B. 

Note that r{x + h) < r{x + h) for each h E Xca^ such that \\h\\ = \\h\\ = 70 and 
x + h€ Int(7rB). Thus (17^ . 117^ and the fact that V{7r^+h n B) = c„^r{x + ft.)™ 
(where c„i is the volume of the unit ball in R"') implies H7.17|l . □ 

Remark. The inequalities in H7.26|l imply that for each x £ d{TTB) and for h 
defined in the proof of Theorem 17. II the following estimate holds 

P{x + h) > c||/i||^. 

This inequality means that if m — 1, P{x) is not difFerentiable at x on the boundary 
d{TrB). 
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7.2. Smoothness of the density P{x) in the interior Int(7ri3). We first refor- 
mulate the definition of P{x) in Ht).7|l . By Ht).3|) and (|6.4|) . for each x e X„„^, we 
see that 

TTa; = © X + TTO, with TTO = {U ® {-au), U e X^} ^ ® Xaa^ ■ 

Therefore the ball 

TTx C] B — {w — u (B V £ TTx : \\w — t&lp < r^}, 

where the center w = w{x) and the radius r = r(a;) are defined in l|7.1(l - (|7.2|l . can 
be rewritten as 

TT^nB = OQx + Vix). 

Here 

V{x) ={w = u® {-au) e TTo : \\w - w{x)\\^ < r^{x)} (7.27) 
with r{x) being defined in H7.2|l and 

w{x) = [{E + a*a)-^a*x] [~a{E + a*a)-'^a*x]. (7.28) 
We now decompose w and w in the basis 9j of ttq, defined in (|6.13|l : 

m m 

and consider the integral 

Pix)= / r(i(;,x)du;= / T{w,x)dw, (7.29) 

Jbdtt^ Jv(x) 

with w = (wi, • • • , Wm) e K™- 

Let a: G Int(7ri?). To investigate the smoothness of P{x) at a;, we consider the 
difference 

P{x + h) - P{x) = I T{w,x + h)dw- i T{w,x)dw. (7.30) 

JV{x+h) Jv{x) 

By (|t).25|l . the function T{w,x) is infinitely differentiable in x and in w. By Taylor 
expansion 

r{w, x + h)= T{w, x) + {T'^{w, x),h) + o{w, X, h), (7.31) 

where 

\o{w,x,h)\<c\\hf (7.32) 

with c depending on {w,x). As usual, in this case, we denote o(w,x,h) as 0(||/i|p). 
Inserting H7.31() into H7.30|l . we conclude that 

P{x + h) — P{x)= I T{w,x)dw — I T{w,x)dw 

Jv(x+h)\V(x) Jv{x)\V(x+h) 

{r'^{w,x),h)dw + o{\\hf). (7.33) 

IVix+h) 

We now calculate the Gateaux derivative of P(x) at x. Denote e — h/\\h\\ and 
I = \\h\\. Divide (|7.33|l by I and take limit as Z — > 0. For the third term in the right 
hand side, we have 

limy/ {r'^{w,x),le)dw ^ {r'^{w,x),e)dw. (7.34) 

1^0 I Jx>[x+le) Jv{x) 
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Figure 1. 



To find the similar limit for the first and second terms in the right hand side, we 
introduce a kind of "polar" coordinates in the sets 'D{x+h)\'D{x) and 'D{x)\'D{x+h) 
when they are not empty. Suppose that 

V{x + h)\V{x) (7.35) 

Let b be the running point on the part of the sphere dB{x + h) which is the part 
of boundary for the set (|7.35|l . Define 

a = w{x)bndB{x), (7.36) 

where 'w{x)h is the vector with end points w{x) and b. Denote by "0 the magnitude 
of the angle Z.hw{x)w{x+h). By the cosine theorem in the triangle l\hw{x)w{x+h)^ 
we see that 

\\w{x + h)b\\^ = \\w{x)h\\^ + \\w{x)w{x + h)\\'^ - 2\\w{x)b\\\\w{x)w{x + h)\\ cos^/j. 

(7.37) 

Define \\'w{x)w{x + h)\\ = p{h) and ||?x'(a;)fe|| = z. Note that || ii;(a; + /i)&|| = 
r{x + h). Then (|7.37|l may be reformulated as a quadratic equation with respect to 
z: 

z^ - 2pz cos i' + p"^ - (x + h) = 

and therefore z — p cos ip + \J r'^{x + h) — p^ sw? ip. We choose the positive sign 
because for ^/j = 0, we should get z = p + r{x + h). As a result, 

||a6('i/;)|| ~ z ~ r{x) — pcosip + \J r'^{x + h) — p^ sin^ ip — r{x). (7.38) 
We introduce polar coordinates in the set 17.35|l : 

V{x + h)\ V{x) 3w^{i^, uj, 7) 

with (jp, u) the spherical coordinate: is the angle in Figure 1 and uj is complement 
spherical coordinate; the coordinate 7 is the distance from a to the point that has 
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coordinate w where a is the fohowing point: a — w{x)wC]dB{x). As is weU-known, 
the Jacobian of the transformation w — w{iIj,uj,j) is equal to {r(x) + 7)™^"'^. So 

dw = (r(x) + -f)""'^d^jdujdj. (7.39) 

Let us estimate ((735|l for smah \\h\\. By (|7!^ . 

p^{h) = \\{E + a*a)-^a*h\\^ + \\a{E + a*a)-'^a*hf (7.40) 

and using (|7.2|l . 

r2(x + /i) = r2(a;) - 2((i; + a*a)"^a*a;, (E + a*a)-^a*h) 

- 2{x - a{E + a*a)-'^a*x, h - a{E + a*a)-^a*h) + 0^(/i). 

Hence by Taylor expansion and by (|7.11() . 

\\abNM{x,h) = p(h)cosij~ -^((E - a(E + a*a)-^a*)x,h) + o(\\h\\^). (7.41) 

r{x) 

Let di{x, h) be the part of the boundary for the set 'D{x + h)\ ^{x), composed 
of the points on sphere dT>{x). Changing to polar coordinates, and applying the 
Taylor expansion for F and using H7.39|l and 17.41(1 . we get 

. A\abm\(x,h) 

T{w,x)dw = / / T{ip,uj,j,x){r{x)+j)"^~^d'ydipdLj 



V{x+h)\D{x) Jdi{x,h)J0 

r(i^, 0, a;)^i(x, V, e)d^dcj||/i|l + 0(||/if ), 

di (x,h) 

(7.42) 

where, recall, e — h/\\h\\, and the function 'I'(x, i/', e) is defined by 

p\\7lM\\{xJi) , _ 

/ ir{x) + iT-^d-i = -{{r{x) + ||a6(V^)|| (x, /i))" - r(x)™) 

io 

= m{x,^,e)\\h\\+o{\\hf). 

In other words, by H7.41|l . 

-^{x.i^.e) = r™-i(a;) (/9(e) cos V'- -^{{E - a{E + a* a)-'^a*)x,e)). (7.43) 

r[x) 

To calculate the integral over 'D{x)\'D{x + h) in (|7.33() we use notations of points 
'w{x), w{x+h), a' , b' and angle ip from Figure 2. Similarly to the case 'D(x + h)\'D{x) 
we define \\w{x)w{x + h)\\ = p{h), \\'w{x)h'\\ — z, \\w{x + h)b')\\ — r{x + h). Then 
by Cosine Theorem in triangle Ab'w{x)w{x + h) we have the equality 



2zp cos tp — r'^{x + h) =0 and therefore z = p cos + \Jr'^{x + h) — p^ sir? ^ (We 
put sign "plus" before square root because for ip — n the equality z = r(x + h) — p 
holds). Then 

ll^ll = r(a;) - z = -p(/i) cos V + r~'^{x){{E - a{E + a*ay'^a*)x, h) 
and similarly to (|7^ . im^ we get 

-/ T{w,x)dw^i V{iP,uj,Q,x)^{x,i),e)diljdw\\h\\+o{\\h\\'^), 

Jv{x)\D(x+h) Jd^{x,h) 

(7.44) 

where ^'(a::, ■0, e) is defined in (|7.43|) and d2{x, h) is the part of the boundary for the 
set V{x) \ V{x + h) composed of the points on sphere d'D{x). 
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Substituting (|7^ . ifT^ . (|7^ into (|73^ . dividing the obtained equation by 
I = ||/i|| and taking the limit Z j 0, we conclude 

^.^ P(x^Ae)^^P(x) ^ /■ r(V',w,0,x)*(x,i^,e)dV'dcc;+ / (r;(w, x), e)dw, 
^1-0 Jav(x) Jv{x) 

This equality shows that the density P{x) possesses the first variation ^ for each x 
on Int(7ri?), and moreover 

P\x,h)^ I T{iP,uj,0,x)'^{x,ip,h)dipdLu + {T'^{w, x), h)dw. (7.45) 

JdV{x) Jv{x) 

Theorem 7.2. The first variation P'{x,h) is the Lagrange variation. Moreover, 
the function 

x^ sup |P'(x, e)| is continuous for X ^ Int{'K B) (7.46) 

Proof. Note that the second terms in right sides of l|7.45|l . (17.43(1 are linear with 
respect to h (to e). Calculation of ((7.43|l in the case of — e gives that the first term 
in right side of H7.43|l is equal r'"^^ (a;)p(— e) cos-i/'i where Vi = ^ + it and ij} is the 
angle from H7.43|l . Hence, 

r™"-^(a;)p(-e) costAi = r™"-^(x)p(e) cos (?/; + tt) = -r'""-^(a;)p(e) cos V'. 

Therefore P'(x, -h) ^ -P'{x, h). 

The assertion (|7.46|) follows directly from the explicite formulas (|7.45|) . H7.43|I . □ 

Now we are in position to prove that P{x) satisfy H3.18|l . 

^Recall (see, for instance, 0|) that P{x) possesses the first variation at a point x if for each 
h g j'fo-CTfc there exists a limit lim\io{P{x + Xh) — P(x))/\ := P'{x,h). Evidently, P'{x,h) is 
positively homogeneous on h : P'{x,Xh) = XP'{x,h), VA > 0. The first variation P'(x,h) called 
Lagrange variation if P'{x, —h) = —P{x, h). 
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Lemma 7.3. For every vi,V2 G ^aak 

j \P{x-Vi)-P{x-V2)\dxi:c\\v^^V2\\x^^^ (7.47) 

where c > does not depend on Wi,W2- 

Proof. For each x, vi,V2 € Xaa^ 

PI \ P( \ r dP{x + V2 + B{v2~vx)) 

P(x - vi) - P[x ~ V2) = J ^ de (7.48) 

Note that derivative dP/dO is weU defined for every 9 £ [0, 1] except, maybe, one 
value such that x + V2 + d{v2 — wi) G dnB because P — outside nB and P 
possesses Lagrange variation inside ttB. Besides, by Theorems 17.21 and 17.11 the 
function 6 dP/dO is integrable and, hence, H7.48|l is well defined. Therefore 
(fnSjl implies 

/ \P{x~vi)- P{x-V2)\dx 

^\\V2-Vi\\x,^^ f f \P'(x + V2+9iv2-Vi), ""^'l^ )\dedx 

"Jx^^^Jo \\V2-Vi\\x,,^ 

^\\v2-vi\\x,,l sup \P\x,e)\dx 

This implies inequality (|7.47|) □ 

Thus, we have proved that RDS (|2.^U|I is a particular case of RDS l|^H5|l and 
therefore it satisfies all conditions of Theorem 13.21 Indeed, Lemma 17.31 implies 
that random dynamical system (|2.31l) satisfies the condition (|3.18l) . As was shown 
above, conditions (|3.16I) . H3.17|I are true for RDS H2.31|l in virtue of Lemma Wa\ So 
assertion of Theorem 13.21 is true for RDS (|2.31l) . This proves Theorem 13.11 and 
completes our investigation in this paper. 
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